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Abstract 

Context. Resonant absorption cyclotron features are a key diagnostic tool to directly measure the strength of the magnetic field of 
accreting neutron stars. However, typical values for cyclotron features lie in the high-energy part of the spectrum between 20 keV and 
50 keV, where detection is often damped by the low statistics from single pointed observations. 

Aims. We show that long-term monitoring campaign performed with Swift-BAT of persistently, but faint, accreting high-mass X-ray 
binaries is able to reveal in their spectra the presence of cyclotron features. 

Methods. We extracted the average Swift-BAT 15-150 keV spectrum from the 54 months long Swift-BAT survey of the high-mass X- 
ray source IGR J16493-4348. To constrain the broadband spectrum we used soft X-ray spectra from Swift-XRT and Suzaku pointed 
observations. 

Results. We model the spectra using a set of phenomenological models usually adopted to describe the energy spectrum of accreting 
high-mass X-ray binaries; irrespective of the models we used, we found significant improvements in the spectral fits adding to the 
models a broad (10 keV width) absorption feature, with best-fitting energy estimate between 30 and 33 keV, that we interpret as 
evidence for a resonant cyclotron absorption feature. We also discuss instrumental issues related to the use of Swift-BAT for this kind 
of studies and the statistical method to weight the confidence level of this detection. Correcting for the gravitational redshift of a 1 .4 
Mq neutron star, the inferred surface magnetic field is B sllrf ~ 3.7 x 10 12 Gauss. The spectral parameters of IGR J16493-4348 fit well 
with empirical correlations observed when the whole sample of high-mass binaries with detected cyclotron features is considered. 

Key words. X-rays: binaries - X-rays: individual: IGR J16493-4348. 



1. Introduction 



O \ The Burst A lert Telescope (BAT , |Barthelmy et alj [200$ on 



OO 

o 



X 



board Swift ( Gehrels et al. 2004 ) has been performing a con- 
tinuous monitoring of the sky in the hard X-ray energy range 
(15-150 keV) since November 2004. The telescope, thanks 
to its large field of view (1.4 steradian half coded) and its 
pointing strategy, covers a fraction between 50% and 80% of 
the sky every day. This has allowed the detection of many of 
the new INTEGRAL Hi gh-Mass X-ray Binaries (HMXBs, e.g. 
|Cusumano et al. 2010b) and the collection of their long term 
light curves and spectra. The long and continuous monitoring of 
these sources allows to investigate the intrinsic emission vari- 
ability, to search for long periodicities (orbital periods) and to 
discover the presence of eclipse events. The role of Swift-BAT 
(hereafter simply BAT) is therefore fundamental to unveil the 
nature and the geometry of these binary systems. Moreover BAT 
collects their long-term averaged energy spectra in the 15.0-150 
keV band. For HMXBs this energy range is extremely important 
as resonant cyclotron scattering features ( CRSFs) are usually ob- 
served in this part of the X-ray spectrum ( Heindl et al. 2004). 

In this paper we analyse the soft and hard X-ray data col- 
lected by BAT, INTEGRAL and by pointed Suzaku and Swift- 
XRT observations on IGR J 16493 -4348 . This source was dis- 



terpart at a posi tion of RA(J2000) = 16°49'26.92"; Dec(J2000) 
= -43°49'8.96"( |Kuiper et al J|2005 ) allowing the association with 
the infrare d counterpart 2MAS S J 1642695-4349090, a B0.5 lb 
supergiant. Mespolietal. (2010) performed K-band spectroscopy 
of the inferred counterpart, confirming the B0. 51 Ia-Ib com- 
panion spectral type; the infrared extinction, translated into the 
equivalent hydrogen column density, was estimated to be (2.92 
± 1.96) x 10 22 cirT 2 , a value that is lower with respect to the 
X-ray absorption, thus indicating that part of the X-ray absorp- 
tion may be local to the compact object. Cusumano et al. (2010a) 
found a periodicity of 6.732 ± 0.002 d, interpreted as the orbital 
period of the system. The folded light curve shows presence of 
an eclipse, lasting ~ 12% of the P or b- Assuming a 1.4 M neu- 
tron star (NS) and 32 R for the companion star, the system has 
mos t likely a low eccentricity value (e < 0.15). Recently, C orbet 
et al. ( pOlCj ) found in RXTE observations evidence for a shorter 
periodicity at 1069 + 7 s, interpreted as the NS spin period. 



A spectral analysis using non simult aneous data from Swift- 
XRT and INTEGRAL ( ftill et al.|p008| ) showed that the broad- 
band X-ray spectrum could be well fitted by an absorbed (Nh 



covered by INTEGRAL in 2004 (Bird et al. 2004) and it was 



initially associated with the radio pulsar PSR J 1649-4349 be- 
cause of a spatial coincidence. A later INTEGRAL observation 
with a deep exposure allowed to reduce the positional uncer - 
tainty and to reject the pulsar association ( Grebenev et al. 2005| ). 
A follow-up observation with Chandra found a soft X-ray coun- 



=5.6 x 10 22 cm" 2 ) p ower-law (F = 0.6± 0.3) with a high-energy 
cut-off at ~ 17 keV. ^/[orris et~aT (2009) analysing Suzaku data 
of IGR J 16493 -4348 , proposed an alternative modelling of the 
X-ray spectrum using a partial covering component multiplied 
by a simple power-law. The covering fraction of the total emis- 
sion was estimated in 0.62±0.07, while the local value of the Nh 
was estimated in ~ 30 x 10 22 cirT 2 . 



1 



A. D'Ai et al.: A cyclotron line in IGR J16493-4348 



2. Data Reduction 

The raw BAT survey data of the first 54 months of the Swift 
mission were retrieved from the HEA S .ARC public archive ^] and 
processed with a dedicated software ( Segreto et al. 201C[ ), that 
performs screening, mosaicking and source detection on BAT 
data and produces spectra and light curves for any given sky po- 
sition. The light curve of IGR J 1 6493-4348 was extracted in the 
15-150 keV energy range with the maximum available time res- 
olution (~ 300 s). The spectrum was obtained by extracting the 
source count rates in 16 energy bands and analysed using the 
BAT spectral redistribution matrix^. The last three energy chan- 
nels (100-150 keV range) were rebinned into a single channel 
in order to have a S/N ratio above 3 <x. The source is detected 
at a significance level of ~ 21 standard deviations. The aver- 
age count rate in the BAT light curve is 1.03 x 10 4 count s _1 
pixel 1 . When considering the light curve at the highest reso- 
lution the maximum of the deviation from the average rate is 
about 7<x, corresponding to an increase in rate of a factor ~ 20. 
Although the count rate indicates brightness variations in the ac- 
creting source, we found no evidence for spectral shape variabil- 
ity, when spectra collected on one month time span were fitted 
using a simple power-law model, with photon indices consistent 
with the average value within the errors. The 15-50 keV and 50- 
150 keV time-averaged fluxes are (3.0+0.3) x 10~ n erg cm 
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s and (I.4+Q3) x 10 erg cm s , respectively. 

We analysed also INTEGRAL-ISGRI (hereafter simply 
ISGRI) long-term data obtained from the Integral high-level 
products archive HEAVENS0. ISGRI data were used in the 
15-150 keV energy range. 

To have the necessary coverage of the soft X-ray band we 
used data of the Swift-XRT (hereafter simply XRT, B urrows 
et al. |2005 ) observation performed on 2006 March 11 (Obs.ID 
00030379002), for a total exposure time of 5.6 ks and a Suzaku 
( Mitsuda et al. 2007 ) obs ervation performed on 2006 October 1 
(Obs.ID 401054010, see |Morris et alj |2009|, for further details). 



The XRT data were processed with standard procedures 
(XRTPIPELINE v.0.12.4), filtering and screening criteria, using 
FTOOLS in the HEASOFT packag e (v 6.8) . The source was ob- 



served in Photon Counting mode (Hill et al. 2004). We adopted 
standard grade filtering 0-12. The source events for spectral 
analysis were extracted from a circular region of 20 pixels ra- 
dius (1 pixel=2.36") centred on the source position as deter- 
mined with XRTCENTROID. The background was extracted from 
an annular region centred on the source with radii of 70 and 130 
pixels, respectively. The spectrum was analysed using ancillary 
response files generated with XRTMKARF and spectral redistri- 
bution matrix vOl 1 . 

We extracted Suzaku scientific data (Obs.ID 401054010) 
from clean events reprocessed with the 2.0.6.13 standard 
pipeline procedure. We obtained four spectra from t he X-ray 
Imaging Spectrometer (XIS(0-4), Koyama et al. 2007| ). Spectra 
from the front-illuminated CCD (XIS0, XIS2 and XIS3) were 
combined into a single spectrum (XIS023) using the addas- 
caspec script, while the XIS 1 spectrum was used in the analysis 
as a different dataset. 

The source extraction region was centred on the CHANDRA 
source position with a radius of 140 ". Background region was 
obtained from an equal spatial region adjacent to the source 



1 http://heasarc.gsfc.nasa. gov/docs/ archive . 
htmT 

2 http : / /heasarc . asfc . nasa . aov^ docs/heasarc/ 
caldb/ data/ swift /bat /index . html 



regio n. Data from the HXD/PIN instrument ( Takahashi et al 
2007 ) were obtained using the hxdpinxbpi script. The HXD/PIN 
background is the sum of a time-variable instrumental back- 
ground (non-X-ray background, NXB) induced by cosmic rays 
and trapped charged particles in the satellite orbit and the intrin- 
sic cosmic background (cosmic X-ray background, CXB). We 
used the 2.0 version of the tuned NXB model sp ectrum released 



by the HXD instrument team (LCFITDT model, Fukazawa et al 



3 http : / / www . isdc . unige . ch/ heavens / 



|2009| ) and the CXB spectrum obtained from the |Boldt| ( |1987| ) 
X-ray background model convolved with the PIN response for 
the flat emission distribution at the epoch of this observation. 
Given the lack of local narrow features, spectra from XIS and 
HXD/PIN were coarsely rebinned into a spectrum of 84 (21) en- 
ergy channels for the XIS023 (XIS1) and 9 energy channels for 
the HXD/PIN. This choice leaves a S/N ratio above 3 for each 
energy channel. XIS data are background dominated below 2 
keV and above 10 keV, so that these channels were not used in 
the spectral analysis. HXD/PIN spectrum was used in the 12.0- 
40.0 keV energy range. 

Spectral fits were performed using XSPEC v.12.6.0. Errors 
are at 90 % confidence level, if not stated otherwise. 



3. Spectral Analysis 

BAT data were collected during the monitoring campaign from 
2004 December to 2009 May, ISGRI data from 2003 January to 
2009 April, and they both track, within the same temporal win- 
dow, the long-term average hard X-ray emission. Swift-XRT and 
Suzaku observations, as indicated in the previous section, are 
pointed observations taken at different times. In order to con- 
strain the broadband spectrum we distinguished, therefore, the 
two soft X-ray observations performed with XRT and Suzaku. 
For Spectrum 1, we hereafter assume the spectrum containing 
the following datasets: XRT, BAT, ISGRI. For Spectrum 2, the 
spectrum composed of these datasets: Suzaku (XIS1, XIS023 
and HXD/PIN), BAT and ISGRI. 

We adopted an inter-calibration multiplicative constant to 
take into account flux uncertainties in the response of the dif- 
ferent instruments. A common fit of the ISGRI and BAT data 
does not reveal any systematic flux mismatch. Data could satis- 
factorily be fitted using a power-law model and the residuals did 
not show any systematic difference in the two datasets. A cali- 
bration constant between the datasets (fixed to one for the BAT 
data and left free for the ISGRI data) was found to be consis- 
tent with one. Therefore, for the following analysis we fixed the 
inter-calibration factor between these two instruments to one. 

For Spectrum 1, the constant was fixed to 1 for the XRT 
dataset and left free for the BAT/ISGRI datasets (C bat )- For 
Spectrum 2, the constant was fixed to 1 for the XIS023 spec- 
trum and left free for the HXD-PIN dataset (CW), for the XIS1 
dataset (C xis i) and for the BAT/ISGRI dataset (C bat )- The 1- 
10 keV absorbed/unabsorbed flux was (5.26+0.08/8.15+0.15) 
x 10~ n erg cirT 2 s _1 for the Swift-XRT observation, whereas 
it was (1.36+0.01/2.35+0.05) x 10" 11 erg cirT 2 s" 1 during the 
Suzaku observation (errorbars take also into account the uncer- 
tainty associated with different modelizations of the spectra). 

To model the data, we employed the three most commonly 
used models to fit broadband spectra of a ccreting high-mass X- 
ray binaries (see ^/[akishima et al. 1999, for a comprehensive 
discussion of the models): a cut-off power law model (model 
cutoffpl), a power-law with a Fermi-Dirac cut-off (model 
f dco), the negative-positive power-law model (model npex, 
with the positive photon index fixed to +2). The cutoffpl 
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model has three main parameters: the photon-index (F), the cut- 
off energy (E cut ), the normalization of the power-law (N po ); the 
f dco model has one more parameter, the folding energy of the 
power-law (Ef i); the npex is basically the sum of two power- 
laws with a common exponential cut-off. 

We did not use the parti a l cove ring model of a power-law 
spectrum as in Morris et al. ( 2009] ), because this model is not 
physically consistent with the nature of IGR J 16493 -4348 , as 
the 6.732 d orbital period and the 1069 s spin period place un- 
ambiguously IGR J16493-4348 in the wind-fed HMXB zone of 
the Pm-b-Pspin diagram (also known as the Corbet diagram, Corbet 



Table 1. Spectral fitting results 



1984). 



3.1. Spectrum 1 

In the first column of Table)!] we present the best-fitting results 



adopting our set of spectral models for the Spectrum 1 datasets. 

The absorbed power-law with a high energy exponential 
cut-off gives a column density Nh = (7.0+1.2) xlO 22 cm" 2 , 
r = 0.7 ± 0.3 and a c ut-off energy jW — 20 + 5 keV. This model 



already proposed in Hill et al. ( 2008 ) is a better representation 
of the data with respect to the more sophisticated f dco model, 
but we found that the npex model provides still a better fit to 
the data, even if the steepness of the power-law, the cut-off en- 
ergy and the Nh are consistent in the two models. Irrespective 
of the model adopted, broad residuals were evident between 30 
and 40 keV (see Fig.]!]). Inspecting the data, we noted that the 
residuals were mostly driven by the residuals of the BAT data, 
because of the better S/N ratio of these data in the hard X-ray 
band. We have tentatively identified the shape of the residuals 
with a broad feature in absorption and we repeated the fits in- 
troducing a cyclotron absorption feature (eye labs component, 
second column in Table [j]), which is expressed according to the 
formula: 



CYCLABS(E) = D c 



(W C E/E C ) 2 
(£-£ c ) 2 + W 2 



(1) 



where E c , D c and W c are the cyclotr on energy, depth and width 
respectively ( Makishima et al. 1990 ). We assume that the broad 
feature constitute the fundamental and, given the low S/N at 
higher energies, we neglect higher harmonics. 

The line parameters could not be all adequately constrained. 
We found that the width of the line strongly correlated with the 
cut-off energy of the continuum models, resulting in a large 
uncertainty. We, therefore, kept this parameter frozen to the 
value of 10 keV, as this was the most reasonable value that pro- 
vided the lowest \ 2 m a ll the models adopted. This choice will 
also be discussed a posteriori. The energy of the line, for the 
npex best-fitting model is 33+4 keV, while the depth of the line 
is 0.5+0.2. Similar values were found also adopting the other 
two models, indicating a weak correlation with the broadband 
shape of the continuum. The addition of this component signifi- 
cantly improved the fit (reduced^ 2 1.16 against a value of 1.32 
without the line). To fit CRSFs, a Gaussi an absorption mode l 
(gabs in Xspec) has been often also used (Coburn et al. 2002| ). 
Replacing the eye labs with the gabs model, did not signifi- 
cantly change the conclusions of our analysis. We found again a 
tight correlation between depth and cr of the line. The best-fitting 
position of the line varied, according to the model adopted, be- 
tween 3 1 keV and 35 keV. Keeping the line <x frozen at 5 keV, 
the best-fitting depth of the line was found between 1 and 5. 
Relative errors on the best-fitting parameters were of the same 
order of the corresponing eye labs analogues. The^ 2 values 
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were found marginally worse with respect to the models with 
the cyclabs component. 

The best-fitting value of Cbat is ~0.2. We verified that this 
value is consistent with the intensity of the BAT light curve dur- 
ing the XRT observation, a factor ~ 4 higher than the average 
source intensity along the 54 months of monitoring. 

We show in Fig.]! data of Spectrum 1 and residuals for the 
models of first/second columns of Table]!]. 

3.2. Spectrum 2 

We used the pointed Suzaku observation to perform a comple- 
mentary analysis, checking how a different soft X-ray spectrum 
could affect the shape of the hard X-ray spectrum. The analysis 
was performed in analogy with the steps described for Spectrum 
1. The third and fourth column of Table]!] show the results of the 
spectral fitting for the three models, with and without the addi- 
tion of the cyclotron line. The results are well in agreement with 
those obtained for Spectrum 1. Again, the best-fitting model of 
the continuum emission is found for the npex model. For this 
model, the normalisation constant Cbat is ~1.4. A value that is 
~6 times higher with respect to Spectrum 7; the observed flux 



3 



A. D'Ai et al.: A cyclotron line in IGR J16493-4348 



Spectrum 1 - DATA 
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Figure 1. Upper panel: datasets that compose the Spectrum 1 
analysis. Black data: XRT spectrum; red data: BAT spectrum; 
green data: ISGRI spectrum. Lower panels: residuals in units of 
cr for the continuum models of Table [l|; upper mid-panels: only 
continuum, lower mid-panels: continuum with the addition of a 
cyclotron absorption feature at ~ 3 1 keV. 
4 



during the Suzaku observation is, in fact, ~4 times less with re- 
spect to the XRT observation. Despite the drop in the soft X-ray 
flux no statistically significant variation in the spectral shape of 
the emission is observed. Best-fitting values for the two spec- 
tra are in agreement, although the power-law steepness is only 
marginally compatible, indicating during the Suzaku observa- 
tion a possible softer X-ray spectrum. 

The best fit to the continuum emission (npex model) gave a 
reduced x 2 = 1-16 (132 d.o.f), while the model with the cyclabs 
component gave a reduced x 2 = 1-10 (130 d.o.f.). The line shape 
is found to depend only marginally on the soft X-ray band. The 
cyclotron line parameters are, however, less constrained with re- 
spect to Spectrum 1, although best-fitting values remain well 
consistent. We show in Fig.B data of Spectrum 2 and residuals 
for the models of third and fourth column of Table [j]. 

3.3. Calibration uncertainties of the BAT instrument 

Because the detection of the cyclotron features is mostly pivoted 
by the BAT high-energy data, we checked for possible biases in 
the analysis due to incorrect background subtraction or system- 
atics in the BAT response matrixn by comparing the extracted 
spectrum of IGR J 16493 -4348 with the spectra of nearby X-ray 
sources having similar spectral shape. 

We selected from the 54-months BAT catalogue the three 
HMXBs closest to IGR J 16493 -4348 : IGR J16393-4643, IGR 
J 164 18-4532 and AX J 1700.2-4220. We extracted the BAT long- 
term spectra and used pointed Swift-XRT and XMM-Newton 
(Epic-PN spectrum for AX J 1700.2-4220) to obtain their broad- 
band X-ray spectra. The spectra of IGR J1639. 1-4641 and AX 
J1700. 2-4220 could be well fitted using a cut-off power-law 
model, while the spectrum of IGR J 1641 8-4532 required a npex 
model to be satisfactorily fitted. In each case the BAT data did 
not show any structure in the residuals in the 30-40 keV range. 
Because of the similarity in the shape of the continuum emission, 
the common background estimate, and the comparable statis- 
tics in the BAT data among all the sources of this sample and 
IGR J 16493 -4348 , we conclude that the residuals shown in the 
spectrum of IGR J 16493 -4348 cannot be due to calibration un- 
certainties of the BAT instrument nor to an incorrect estimate, or 
subtraction, of the X-ray background. 

3.4. A method for testing the significance of the detection 

To estimate the statistical significance of the residuals present in 
the BAT data of IGR J 16493 -4348, we adopted the following 



procedure (see also Suchy et al. 201 1 for a similar approach) : 



we assumed the best-fitting npex model for Spectrum 1 
in Table [j] without the cyclotron line as the null-hypothesis 
model; 

we simulated accordingly a faked BAT spectrum having the 
same S/N ratio of the real BAT spectrum; 
we fitted again the datasets of Spectrum 1 (real data and sim- 
ulated BAT spectrum) using the npex model with and with- 
out the addition of an absorption cyclotron line (width con- 
strained at 10 keV, with line energy and optical depth as free 
parameters); 

we calculated the difference in A^ 2 for the two best-fitting 
models: npex and npex+cyclabs; 



4 http : / /swift . crsf c . nasa ■ gov/docs/ swift/ 
analysis/bat_digest . html 
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Spectrum 2 - DATA 
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Figure 2. Upper panel: datasets that compose the Spectrum 2 
analysis. Black data: XIS023 spectrum; red data: BAT spec- 
trum; green data: ISGRI spectrum; blue data: Suzaku/HXD-PIN 
spectrum; light blue: XIS0 spectrum. Lower panels: residuals 
in units of cr for the continuum models of Table [|; upper mid- 
panels: only continuum, lower mid-panels: continuum with the 
addition of a cyclotron absorption feature at ~ 3 1 keV. 



Figure 3. Histograms of the difference in Ax 2 obtained fitting 
the datasets of Spectrum 1 and Spectrum 2, when the real BAT 
spectrum is substituted with a BAT simulated spectrum. The x 2 
difference is calculated for the two best-fitting models: npex and 
npex+cyc. 
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- we repeated this procedure for other 10,000 faked BAT spec- 
tra; 

- we plotted with histograms the Ax 2 difference versus the 
number of spectra; 

- we repeated the same procedure for the Spectrum 2 datasets; 

In FigJ|, we show the results of our method. In 10,000 spec- 
tra we did not obtain any difference in the Ax 2 as the ones re- 
ported in Table [j], both for the Spectrum 1 (observed A^f 2 =19) 
and for the Spectrum 2 (observed A^ ,2 =10) datasets. The shape 
of the distribution is different in the two cases: the XRT data 
constrain less the broadband shape and the addition of a broad 
feature in absorption in most cases improves the fit. On the con- 
trary the Suzaku data are less sensible to fitting improvements 
even for the addition of a broad, continuum-like, feature. This 
is expected since the statistical weight of the Suzaku XIS (~ 
20,000 counts in the XIS spectra versus the ~ 2,000 counts of the 
XRT spectrum) data more strongly pivot the overall determina- 
tion of the continuum shape and the x 2 result. At the same time, 
tighter constraints of the continuum shape make random fluctu- 
ations as artifacts less probable. We note, therefore, that the de- 
tection level as estimated from the difference in the^ 2 of the two 
models (with and without the addition of the cyclotron absorp- 
tion feature) in Spectrum 2 real data, although in absolute value 
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less than what found for Spectrum 1, is consistent and in agree- 
ment with the detection level obtained from Spectrum 1 data. 
Moreover, we have estimated this detection level for the model 
(npex) that provided the lowest improvements for the addition 
of this features; the other spectral models would have given a 
much higher detection confidence. Because from the generation 
of 10,000 simulated spectra no random fit improvement of the 
same order, or greater, than the one we obtained from real data 
was obtained for both datasets, we estimate a detection level cor- 
responding to at least ~3.9 gaussian sigmas for the presence of a 
cyclotron line in our fits. 

4. Discussion 

There are 15 sources that clearly have shown CRSFs signatures 
in their spectra. This number represents about 10% of the total 



population of HMXBs detected in our Galaxy ( |Reig|[201 1[ ) and, 
despite the considerable amount of observing time dedicated by 
modern and past X-ray observatories, to this class of sources, 
new detections of CRSFs have not considerably risen in the last 
decade. 

We have presented new results based on the analysis of the 
data collected by BAT during the first 54 months of the Swift 
mission on a persistent sgHMXB IGR J 16493 -4348 . We have 
shown that the long-term BAT collected hard X-ray spectra pro- 
vides an access for detecting the presence of CRSFs in this kind 
of sources. This new approach relies on a series of assumptions 
and caveats that will be hereinafter discussed. 

The broadband spectrum, modelled with a positive-negative 
cut-off power-law model provides the statistically most favoured 
description of the data in the examined datasets. The residuals in 
the BAT energy range show the presence of an absorption feature 
between ~28 and ~38 keV. The feature is not sensibly dependent 
on the adopted soft X-ray contribution, although we have just ex- 
ploited the only two available to date pointed observations with 
XRT and Suzaku. The broad negative residuals are compatible 
with a resonant cyclotron absorption feature with an estimated 
significance of the detection at least at ~ 4 cr level. 

Similar features are often seen in the high energy spectra of 
sgHMXB. They are interpreted as cyclotron lines produced near 
the magnetic poles of the accretion-powered neutron star. These 
are due to resonant scattering processes of the X-rays by elec- 
trons whose kinetic energies are quantised in discrete Landau en- 
ergy levels perpendicular to the B-field. Their detection is funda- 
mental for the understanding of the accretion mechanisms unto 
the neutron star as it allows a direct measurement of the magnetic 
field of the neutron star (or at least a lower limit, depending on 
the height where the resona nt line absorption is effectively pro- 
duced, 



see 



NishimuK 2011) 



We consider here that the cyclotron absorption takes place at 
the poles of the neutron star, taking into account the relativistic 
gravitational redshift: = E cyc (l + z) , with 



(\+zT l = 1- 



2GM X 
Rxc 2 



0.5 



(2) 



Assuming typical values for the NS mass and radius (1.4 
M and 10 km, respectively) and the observed value of the cy- 
clotron line energy of 33+4 keV (npex model for Spectrum 1), 
this implies a magnetic field of Bns (3.7+0.4) xlO 12 Gauss at 
the surface of the NS. 

The energy resolution of the spectra is not sufficient to inde- 
pendently constrain line energy, line width and depth. Because 



the line energy and depth could be relatively better constrained 
with respect to the line width, we choose to freeze this value 
to 10 keV. This value represents a first-order guess assumed in 
consideration of the emerging relations among the parameters 
that characterize the cyclotron line shapes. An empirical corre 



lation between line energy and line width (Coburn et al. 2002) 
is, in fact, observed, where cyclotron lines detected between 30 
keV and 40 keV show preferably line widths in the 5-10 keV 
rang e (e.g. in MX 0656- 72, 4U 0352+309, XT E J 1946+274, 



[see |Coburn et al] ( |2002[ ); |McBride et al| ( ^006[ )]). We tried in 



all the fittings also to freeze the line width at 5 keV and 8 keV 
to check how sensible was the determination of this parameter. 
We generally found that the 10 keV guess always provided the 
lowest^ 2 , even if the improvements in most cases were marginal 
(2< A^ 2 <5). Another important empirical correlation concerns 
the cut-off energies of the continuum emission and the cyclotron 
line energy, where a correlation is observed with cut-off energies 
typically at half the value of the cyclotron energy. Also in this 
case, our results appear in agreement with the general trend that 
is observed, with a cut-off energ y at ~15 keV and th e cyclotron 
energy at ~ two times this value (Heindl et al. 2004). 



Possible biases in this analysis are envisaged by the use 
of long-term time averaged spectra for features that may show 
some variability. Some luminosity-dependent shifts of the cy- 
clotron line position were in fact clearly ob s erved in some 
sourc es ( Mihara et al. 2004 ; Staubert et al. 2007 ; [Nakajima et al 



2010). However the large uncertainties in the present analysis 
estimates do not allow to set this possibility into investigation. 
IGR J 16493 -4348 has never shown unusual periods of strong 
changes in its luminosity, nor any outbursting behavior as in the 
case of X0331+53, wh ere luminosity changed by a factor ~ 200 
(Nakajima et al. 2010), so that we argue that possible shifts, if 



present, should be within our quoted error-bars. Moreover, it is 
to be noted that the cyclotron feature at 45 keV of A0535+262 
did not show any variati ons despite a chan ge in luminosity of 
two orders of magnitude (Terada et al. 2006). Another source of 



bias concerns the possibility of continuum spectral changes, that, 
when averaged, could result in spectral artifacts. The analysis of 
the soft X-ray spectrum from the XRT and Suzaku pointed ob- 
servations are not sufficient to test any variability pattern in the 
broadband spectral properties; use of simplified, phenomenolog- 
ical, models are also possible source of bias and a more physi- 
cally based model should be employed to also i nfer the physical 
roperties of the accretion co lumn environment ( Becker & Wolff 
007 ; Schonherr et al. 2007 ) and to better constrain the shape of 
cyclotron resonance features. 



4.1. Conclusions 

We have presented spectral analysis in the high-energy X-ray 
band of IGR J16493-4348 using long-term survey data from 
BAT and ISGRI, complemented with pointed soft X-ray obser- 
vations from XRT and Suzaku. The BAT spectrum is to date the 
highest S/N spectrum of this source in the high-energy X-ray 
band. When the spectrum is complemented with soft X-ray data, 
a broad absorption feature is detected, irrespective of the broad- 
band model used to fit the data. We modeled the feature assum- 
ing that it is a CRSF This interpretations is statistically accept- 
able and the physical interpretation of the parameters is in agree- 
ment with what expected from the whole sample of HMXBs that 
show similar features. We estimated the significance level of the 
detection > 3.9 cr confidence level and next generation of future 
observatories dedicated to the high X-ray band (e.g. NuSTAR 
and Astro-H) will likely provide a more stringent confirm of its 
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presence and a substantial improvement in the determination of 
its shape. 
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